On the Physics of the Riemann Zeros 
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Abstract. We discuss a formal derivation of an integral expression for the Li coefficients 
associated with the Riemann ^-function which, in particular, indicates that their positivity 
criterion is obeyed, whereby entailing the criticality of the non-trivial zeros. We conjecture 
the validity of this and related expressions without the need for the Riemann Hypothesis 
and discuss a physical interpretation of this result within the Hilbert-Polya approach. In 
this context we also outline a relation between string theory and the Riemann Hypothesis. 



1. Introduction: Riemann Zeros and the Li Criterion 

This talk is devoted to the physics of the non-trivial zeros of the Riemann ^-function, the 
so-called Riemann zeros [Ij. First, we review Li's Criterion, with an emphasis on requisite 
properties of the Riemann £- function. In the central section of the talk, we show how using the 
exact formula for the number of non-trivial zeros N(T) up to a given height in the critical 
strip, and its derivative as eigenvalue distributions in the critical strip, would imply the 
positivity of Li's coefficients. Finally, we interpret this result from a physicist's perspective 
emphasizing a possibly exciting connection to string theory [2|. 

Let us, perhaps more for the sake of notation, first remind the reader of the statement of 
the Riemann Hypothesis [3H5]: The analytic continuation, from z £ R>i to the whole complex 

oo 

plane z € C of := Yl n~ z has all its zeros in the critical strip S := {z | < Re(z) < 1} 

71=1 

lying on the critical line £ := {z \ Re(z) = ^}. Now, £(z) obeys a remarkable functional 
equation £(z) = 2 z tt z ~ 1 sin(^7rz)r(l — z)C,{l — z), which inspired Riemann to define the 
^-function, which will be central to our discussions: 

i(z):=\z(z-l)^T( Z -)Q(z) . (1.1) 



There are many advantages to considering £(z) instead of C( z )> which we now summarize. 



1.1. The Riemann ^-function 

First, note that C,{z) has trivial zeros at all negative even integers, which are conveniently 
cancelled by the corresponding simple poles of the T-function. Hence, has only non- 
trivial zeros, all located within the critical strip, and affords an elegant Weierstrafi product 
expansion, known as the Hadamard product: 

pes P 

where £(0) = |, and p are the non-trivial zeros of Q{z) in S since none of the other factors in 
the definition of £(z) vanishes therein. It is an obvious but important fact that S extends to 
both the upper and the lower half-planes and that all the p occur in conjugate pairs above 
and below the real line. We shall denote 5+ as the critical strip above the real line and 
correspondingly £ + as the upper critical line. 

Second, the functional equation becomes particularly simple: 

e(*)=£(i-z) . (i.3) 

Indeed, the symmetry about the critical line of Re(z) = I becomes manifest. Hence, not only 
are the zeros symmetric about the real axis, they are also symmetric about the critical line. 
The Riemann Hypothesis postulates that all the zeros in fact lie on C It is known that the 
first 10 13 zeros lie on the critical line [6]. 

Third, consider the conformal mapping z i— > -zj, which maps £+ to the boundary circle 

\z\ = 1 of the unit disk and the entirety of the left half of <S+, i.e., {z | < Ke(z) < i}, to 
the interior of the open unit disk. Consider, therefore, the function 

0(z) := = £(_L) , (1.4) 

where we have used (jl.3p in the second equality. Since all zeros are symmetric around the 
critical line, to consider the zeros of £ within S therefore amounts to considering the zeros 
of 4>(z) within the open unit disk. Indeed, the Riemann Hypothesis would require that there 
be no such zeros and indeed that all the critical zeros are on the boundary circle. Since the 
number of zeros of an analytic function f(z) in a region R is equal to j R dz 4t^t, we have 



/(*) 

that: The Riemann Hypothesis is equivalent to -^jjj being analytic (i.e., holomorphic and 
without poles) within the open unit disk \z\ < 1. 

Fourth, a beautiful exact result relates the number of zeros up to height T within S and 
the ^-function: 

N(T) := #{p G S , < Im(p) < T\ ((p) = 0} = - Im log ^ + iT) . (1.5) 

Extensive work has been done in estimating (jl.5p . We have, using the definition (jl.ip for the 
first equality, that 

N(T) = U(T) + 1 + - Im log cd + ^) ~ ^ log ^ - ^ - J + 0(log T) , (1.6) 

7T 7T Z Z7T Z7T Z7T O 

where historically 1?(T) = Imlogr(|T + \) - | log 7r is known as the average part and 
^ImlogC(^ + iT) ~ O(logT) is the fluctuating part around the essentially TlogT growth 



of the former. This is the inspiration behind classical (the average) and quantum mechanical 
(the fluctuation) interpretations of the critical zeros [TIE] • 

We note that N(T) is a real step function, increasing by unity each time a new critical 
zero is encountered: 

oo 

N(T)= Y J 9{T--imp) = Y t a j 9{T- l i j ). (1.7) 
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It is vital to explain the above in detail. We can explicitly write the upper critical zeros as 
Pj,k = &j,k + indexed by j, k € Z>o, where Oj t k € (0,1) and pj € K>o- The zeros are 
ordered so that Pj+\ > Pj- Crucially, as we do not assume the Riemann Hypothesis, the real 
part of p need not be i. Furthermore, we do not assume simplicity of the zeros, which is as 
yet also unknown [9]. If, for example, we have a double root, we explicitly count this twice. 
By the functional identity, zeros of the ^-function not on the critical line are paired within 
the upper critical strip: if a + ip is a zero, then so is (1 — a) + ip. The index k = 1, . . . , ay 
enumerates the zeros with the same imaginary part. We order the zeros so that crjk+i > <fj k- 
The ay then counts the number of zeros with imaginary part pj including the multiplicities 
of the roots. If, for example, there are a pair of simple roots with imaginary part pj, then 
in an epsilon interval around pj, the counting function N(T) jumps by two. Contrariwise, if 
p € £+ is the only simple root with imaginary part pj, then N(T) jumps by one. We know 
that the number of roots with imaginary part in the interval (0, T) is finite. Indeed, the 
asymptotics of the expression are given in (|1.6j) . In summary, N(T) is a strictly increasing 
step function as we move up in height T regardless of the Riemann Hypothesis. The Riemann 
Hypothesis is, of course, the statement that Oj = \ for all zeros. 

The rewriting of N(T) lends itself to a wonderful interpretation: N(T) is a cumulative 
density function defined over the critical strip. In other words, its derivative, p(T) := N'(T) 
is a density function of distributions. That is, one could conceive of a physical system whose 
energy levels (eigenvalues of the Hamiltonian) are thus distributed; this is along the school of 
thought of Hilbert-Polya [TJ. Of course, the resulting eigenvalue density is highly non-smooth, 
but is, rather, a sum of delta- functions: 

oo 

p(T):=N'(T)=J2^ j S(T-p j ). (1.8) 



1.2. Li's Criterion 

Another striking property of the ^-function was noted by Li [10] not too long ago. Let {k n }, 
for positive integers n, be defined by 

1 c\ n 

then we have that: The Riemann Hypothesis is equivalent to the condition that k n > for all 
n 6 Z>o- This is Li's Criterion. Li also showed two equivalent ways of writing these numbers, 
namely 



<//(z) 
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Y,k n+ iz n ; (1.10) 

(1.11) 
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where we recall that 4>(z) := j). 



2. Density Function for the Distribution of Critical Zeros 

Thus armed with all necessary ingredients, let us proceed to rewrite the Li coefficients in a 
suggestive form. We will do so formally and point out the subtleties involved at the end of 
the manipulations pQ. 

First, integrating (jl.lOp and using the definition as well as the functional equation (|1.4p . 
we have that, for z in the unit disk, 



i 00 r oc , 

log£(-^r) =log£(-— ) = J> n / d^"- 1 = -log2 + ^^^, (2.1 
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where — log 2 is easily checked to be the constant of integration. 

We note, however, that the above expansion has radius of convergence strictly less than 
unity, whereby excluding the unit circle, to which, crucially, C is mapped under our conformal 
transformation z i— > Therefore, it is imperative to analytically continue. Let us start 

with (|2.ip and rewrite the expansion about the point z = — 1. Note that the region of 
convergence for this is an open of circle of radius two centered about z = —1, which, in 
particular, encloses the entirety of the closed unit disk, except the point z = 1, where there 
is a pole for <p(z) = Therefore, we have that 

1 i 00 t, 00 

log£(— ) =log- + ^-." = 6o + E b ^ z + l ) n > ( 2 - 2 ) 

-L Z Zi Tl 
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where we have written the two expansions together for comparison. 

Expanding (z + l) n , we readily obtain an expression for the Li coefficients in terms of the 
new expansion coefficients b n (studied further in reference [H]), k n = nY^ =n (Jjbj • We 
now wish to solve for the k n explicitly and demonstrate positivity. 

We first see, using the counting formula (jl.5p and the expansion (|2.2p . that 



iVOu) 



-Imlog4+v) = E-Imf^ + iy , (2.3) 

n=l x ' 



where we have used the substitution | +ifi = j^, or z = |^rf • Then, since 



Im (^±1 + l) = „ sin(n tan" 1 — ) = 2™ cos™ 9 sm(n9) , cos := ; ^ , 

V2/i-i y (v + i)f v w y ; ' vvti 

(2-4) 

and due to the definition of the Chebyshev polynomial of the second kind 

?7 m _i(cos0) := — — — , (2.5) 
sm m 



we arrive at the conclusion that for all m G Z 



>o: 



:>,2 " l/ d ^ (v + 1)2 ^^i^)' ( 2 - 6 ) 



This is a beautiful and remarkable integral expression for the Li coefficients which works 
wonderfully numerically pQ. 

Let us now proceed to simplify the integral (|2.6p . First, recall that Chebyshev polynomials 
of the second kind are related to the Chebyshev polynomials of the first kind by the equation 

dx U n (x) = — 5—T n+1 (x) . (2.7) 
n + 1 

Then using the integral expression for k m we have that [lj 

00 4/x 2 - 1 

k m = 2j2<Xj (1 ~ T m ( i )) , meZ >0 . (2-8) 
3=1 ^ + 

Recalling the definition of the Chebyshev polynomial of the first kind, that T m (cos9) := 
cos(m8) and remembering that all atj, regardless of actual value, are strictly positive integers, 

each summand in (|2.8|) is non- negative. In fact, because 4 $ +1 < 1 strictly for any imaginary 

part as we know the first fij is at least 14, 



T m {— ) < 1 for any given m € Z>o and for all j G Z>o , (2-9) 

ifj,j + 1 



so each summand is actually strictly positive. In conclusion, this would imply that all Li's 
coefficients k m > 0. Indeed, we have numerically evaluated [T] the first few k m coefficients 
given in (12. 8j) . using the first 10 4 critical zeros, and see that they are very close to the results 
obtained from (jl.lOp and (|2.6p . and moreover constitute an increasing sequence of strictly 
positive numbers. In particular, the above discussions would force all critical zeros to lie 
entirely on L. 

We seem to have arrived the positivity of the Li's coefficients without the assumption of 
the Riemann Hypothesis! This, of course, is not quite true and the criticality of the zeros 
of the ^-function had in fact been inherently invoked in the choice of a branch cut for the 
logarithm. Nevertheless, we have formally arrived the integral formula (j2.6|) which does not 
seem to depend on the precise location of the zeros and the subsequent choice of branch-cut. 
Note that Li's criterion can be applied to other zeta functions [10] and thus it is natural to 
conjecture that our integral and summation formulas for the Li coefficients should exist in 
those cases as well. 



3. The Physics of the Critical Zeros and String Theory 

Attempts to use a physical argument to explain the Riemann Hypothesis date back to Hilbert 
and Polya Were one to find a matrix 

H=^-I + i-T, (3.1) 

whose eigenvalues are all the zeros of £(z) and if in addition the matrix T were Hermitian, then 
all the zeros would lie on the critical line C Therefore, one would wish to find the appropriate 
operator H and its corresponding eigenfunctions, or, in physical terms, the Hamiltonian and 
its wave functions. 



What can one say about the quantum mechanical picture pQ? Let us investigate the 
secular equation of the operator in (|3.ip . It should by construction have roots corresponding 
to the zeros of that is, 

£0) = det(> • I - H) . (3.2) 

To arrive at the appropriate matrix description, we start by considering the expansion (|2,ip . 
which, combined with (|2.8p . gives 



1 oo oo „ 

ioge( T -^) = -io g 2 + x:E — 

1 — s ^-^ n 

n=l j=l 

4M? -1 



T M 



+ 



V 



(3.3) 



Put cos 9j : — §— t . Using the definition of the Chebyshev polynomial and reversing the 



sums, we obtain 



1 OO OO , v 

lo ^— ) = - lo g 2 +E 2 ^E ( - - - cos K) ) * n - ( 3 - 4 ) 

j=l n=l ^ ' 

Substituting z = and using that e ±l9j = 1 — (| ± ifij)" 1 , we find 

which is simply a refined version of the Hadamard product, a reassuring check. Note that 
each factor in the product obeys the functional identity. 

Equation ()3.5[) is to be compared with (j3.2|) . We see that there exists a matrix such that 
its eigenvalues coincide with which is consistent with the Hilbert-Polya picture. As the 
expression vanishes for zeros on the critical line, the refined Hadamard product is already in 
the form of the determinant of an infinite matrix, and the eigenvalues of T are real. 



3.1. Topological String Theory, Gromov-Witten Invariants and the Critical Zeros 
The profound proofs of the Riemann Hypothesis for the case of finite fields by Weil, Deligne 
and others are essentially geometric, and thus one might expect a deeply geometric resolution 
of the Riemann Hypothesis |12j . Motivated by this expectation in this concluding section we 
present a possibly exciting connection between Gromov-Witten (GW) invariants, topological 
string theory and the Riemann zeros [2]. 

First we recall a reformulation of the Riemann Hypothesis in terms of an elementary 
problem due to Lagarias [13]: Let H n = ^j=i J be the n-the harmonic number and let 
<Ti = Yld\n d * s ^ ne sum °f divisors (Ramanujan) function, then the Riemann Hypothesis is 
true if and only if 

<ri(n) < H n + e Hn logH n (3.6) 

for n > 1. In the following we notice the appearance of o~i{n) in the genus-one free energy of 
a particular topological string theory [2]. 

We start with the following statement from reference [13]: Take the type II string on the 
local Calabi-Yau threefold defined by 



uv = W(y,z) = y 2 -P n {z) 2 + l , 



(3.7) 



where W(y,z) = defines a Riemann surface and P n (z) is a degree n polynomial. The 
genus-one amplitude can be given exactly (due to Seiberg-Witten (SW) theory) by 

Fi(r) = -log77(r) ; (3.8) 

this is also the genus-one free energy of the topological string (i.e. topologically twisted from 
the M = 2 SW theory). From the definition of the Dedekind eta-function, we can expand, 
setting q = exp(27m~), 

_l oo oc oc ij 

F 1 (r) = -l og r,(r) = -^-logq-J2^g(l-q i ) = -^--J2J2- 

1=1 8=1 j=l 3 

OO OO j n OO oo 

-^-EE7 = -^-E(Ev = -^-E^' ( 3 - 9 ) 

j=l j=l n=l d|n n=l 

Thus, the GW invariants at genus g = 1 and number of marked points m = 0, for the above 
Calabi-Yau X is the Ramanujan <j\ function! That is: 

N g=1:m=0 {n;X) =(ri(n) . (3.10) 

Then, the statement, due to Lagarias' theorem [T3], is that the asymptotics of the free 
energy of this particular topological string at genus-one world-sheet embedding should be 
(T\{n) < H n + e Hn logH n , if and only if Riemann Hypothesis holds. Therefore: (1), we 
know that the Lagarias bound on a\{n) in terms of the harmonic numbers is equivalent 
to Riemann Hypothesis and (2), we know of examples of GW invariants which are equal 
to a\{n). Hence the question: does there exist an independent geometric bound on such 
invariants? (At this point one might wonder whether the GW invariants could encode the 
critical zeros. This seems unnatural because GW invariants are integers indexed by integers 
whilst the critical zeros are not, whence it is hard to imagine, though no entirely impossible, 
that a combinatorial problem should be associated thereto.) 

Also, by considering rational curves of degree d on the K3 surface we can deduce that [2] 

k 

T{n)N\{k—n) < fccxi(n), for each k = 0, 1, 2, ... , and thus restate the Riemann Hypothesis 

n=0 

in the following form [2]: 

k 

^2 T(n)Nx(k — n)< k(H k + e Hk logH k ), k = 0,1,2,... (3.11) 

n=0 

where r(n) is the famous Ramanujan multiplicative tau-function (for which, r(0) = 0), and 
Ni(j) is the number of elliptic curves with j nodes and passing through a given point on a 
generic, smooth K3 surface. Note that according to reference |15j there exists a universal 
constant a such that 

Rat(X) < X (X) + a , (3.12) 

for every generic smooth surface X, with x(-X') being the Euler number. We need a version 
of this statement for elliptic curves. 

The papers on the so-called Calabi-Yau differential equations have studied the integer 
coefficients of mirror maps, which we know to be an interesting phenomenon in and of itself. 



(In particular, the mirror symmetry might be related to the Birch and Swinnerton-Dyer 
conjecture [IB].) Reference [IT] is particularly interesting. We know that, after judicious 
transformation, the coefficients of the mirror map are precisely the GW invariants. In 
reference [17J, an expression of the coefficients in terms of harmonic numbers is given. The 
strategy is clear. There appear to exist two directions: one giving the GW invariants in terms 
of ci(ra) and the other, the GW invariants via the mirror map, in terms of the harmonic 
numbers. One thus must relate the two and the question is whether the relation is strong 
enough to give the Lagarias bound. 
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